NASA TECHNICAL 

MEMORAND U M 



NASA TM X-2646 









f .Hi: :4 m 







mm? 


COMPARISON pi NEUTRON SPECTRA 

measured write threeIsizes 1 

OF ORGANIC LIQUID SCINTILLATORS 
USING DIFFERENTIATION ANALYSIS 

dt. ' k > >•- i s 


by Donald E Shook and Clarence R. Pierce 

Lewis Research Center 
Cleveland, Ohio 44B5 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • NOVEMBER 1972 


1. Report No. 

NASA TM X-2646 

2. Government Accession No. 

3. Recipient's Catalog No. 

4. Title and Subtitle COMPARISON OF NEUTRON SPECTRA MEASURED 
WITH THREE SIZES OF ORGANIC LIQUID SCINTILLATORS 
USING DIFFERENTIATION ANALYSIS 

5. Report Date 

November 1972 

6. Performing Organization Code 

• 

7. Author(s) 

Donald F. Shook and Clarence R. Pierce 

8. Performing Organization Report No. 

E-6883 



10. Work Unit No. 

9. Performing Organization Name and Address 


503-10 

Lewis Research Center 

National Aeronautics and Space Administration 

11. Contract or Grant No. 

Cleveland, Ohio 44135 


13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 


Technical Memorandum 

National Aeronautics and Space Administration 
Washington, D. ... 

14. Sponsoring Agency Code 

15. Supplementary Notes 

16. Abstract 


Proton recoil distributions were obtained by using three organic liquid scintillators of different 
size. The measured distributions are converted to neutron spectra by differentiation analysis 
for comparison to the unfolded spectra of the largest scintillator. The approximations in- 
volved in the differentiation analysis are indicated to have small effects, on the precision of 
neutron spectra measured with the smaller scintillators but introduce significant error for the 
largest scintillator. In the case of the smallest cylindrical scintillator, nominally 1. 2 by 
1. 3 cm, the efficiency is shown to be insensitive to multiple scattering and to the angular dis- 
tribution of the incident flux. These characteristics of the smaller scintillator make possible 
its use to measure scalar flux spectra within media where high efficiency is not required. 


17. Key Words (Suggested by Author(s)) 

Neutron spectra; Organic scintillator; Neutron 
flux measurement; Spectra unfolding; In-core 
spectrum measurement 

18. Distribution Statement 

Unclassified - unlimited 

19. Security Classif. (of this report) 
Unclassified 

20. Security Classif. (of this page) 
Unclassified 

21. No. of Pages 
28 

22. Price* 

$3.00 


* For sale by the National Technical Information Service, Springfield, Virginia 22151 























COMPARISON OF NEUTRON SPECTRA MEASURED WITH THREE SIZES 
OF ORGANIC LIQUID SCINTILLATORS USING 
DIFFERENTIATION ANALYSIS 
by Donald F. Shook and Clarence R. Pierce 
Lewis Research Center 

SUMMARY 

Proton recoil pulse height distributions obtained by using three organic liquid scin- 
tillators of different size are converted to neutron spectra by differentiation analysis. 

The approximations involved in this analysis are indicated to have small effects on the 
precision of neutron spectra measured with the smaller scintillators but introduce sig- 
nificant error for the largest scintillator. The data were obtained by using three cylin- 
drical scintillators nominally 1. 2 by 1. 3, 2. 2 by 2. 4, and 5. 0 by 5. 0 centimeters. 

For the largest scintillator, an unfolding analysis employing calibrated mono- 
energetic response functions was used as a reference. All four measurements agree on 
the measured shape of a polonium -beryllium neutron source spectrum, but absolute 
fluxes obtained with the largest scintillator by differentiation analysis are 10 to 20 per- 
cent lower than the fluxes obtained by using the calibrated unfolding analysis. However, 
absolute fluxes obtained with the two smaller scintillators by differentiation analysis 
agree with fluxes measured by using the calibrated scintillator. This agreement indi- 
cates that for small scintillators errors in the analytic expression for efficiency and 
multiple scattering corrections used in the differentiation analysis are insignificant. 

In the case of the smallest scintillator, the efficiency is shown to be insensitive to 
multiple scattering and to the angular distribution of the incident flux. This character- 
istic of the smaller scintillator makes possible its use for measuring scalar flux spectra 
within media where high efficiency is not required. The resolution of the smallest scin- 
tillator is indicated to be 50 percent better than that of the largest. The computer code 
TUNS, which was written to analyze proton recoil pulse height data by the differentiation 
method, is described. 



INTRODUCTION 


The method of measuring fast neutron spectra by analysis of proton recoil pulse 
height distributions obtained with organic scintillators has received considerable atten- 
tion (refs. 1 and 2). Using organic liquid scintillators and new electronic circuits 
which permit separation of gamma ray signals allows a differential proton recoil spec- 
trum to be obtained which shows significant spectral detail. A number of methods are 
available for converting the measured proton recoil spectrum into the incident neutron 
energy spectrum, and it is the differentiation method of analysis that is primarily con- 
sidered in this report. This method was first applied in the early 1960’s (ref. 1), but 
was not widely used because (1) instrumentation in use at that time did not permit 
measurement of highly accurate proton recoil spectra, (2) the differentiation method 
involves a calculation of scintillator efficiency with corrections for multiple scattering 
and wall effects that are difficult to evaluate, (3) the method was not verified by direct 
comparison with precise calculations or other methods of measurement, and (4) an ex- 
perimental problem due to the light output anisotropy of stilbene crystals existed. 

More recently proton recoil spectra have been measured by using liquid scintilla- 
tors in conjunction with monoenergetic response functions and the proton recoil data un- 
folded (ref. 2). The method of differentiation analysis is reexamined in this report in 
the light of the recent developments because it offers more flexibility than does unfold- 
ing analysis. 

In this report proton recoil distributions are obtained by using three sizes of organic 
liquid scintillators. The consistent set of absolute measurements of the same polonium- 
beryllium (PoBe) neutron source with different sized scintillators was taken to test the 
differentiation method of data reduction in which the approximations used are expected 
to be dependent on scintillator size. For reference, a comparison is also made with an 
unfolding analysis for the measurements made with the 5- centimeter-diameter by 
5-centimeter-high scintillator for which calibrated monoenergetic response functions 
are available (ref. 2). The unfolding analysis is not subject to uncertainties due to 
multiple scattering and wall effect corrections since a large set of measured mono- 
energetic response functions for the specific scintillator are used. If the spectrometer 
is employed in exactly the same geometry as the response functions have been meas- 
ured, the unfolded spectra should be accurate, although the error in the unfolded spectra 
due to the error in the response functions is difficult to evaluate. The accuracy of the 
unfolding method has been tested by comparing neutron spectra by both the unfolding and 
time-of-height techniques with reasonable success (ref. 3). The characteristics of the 
smallest scintillator used are emphasized in this report. In particular, the improved 
resolution and applications of this scintillator, as summarized in reference 4, are 
discussed. 
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EXPERIMENT 


The measurements were made by using a PoBe neutron source circulated by the 
Oak Ridge National Laboratory. It was positioned 1 meter from the scintillator and 
about 2 meters from the floor at the approximate center of a 4- by 10- by 12 -meter 
room. The three organic liquid scintillators used were contained in glass cylinders 
with the following inside diameters and lengths (cm): 1. 22 by 1. 27, 2. 16 by 2. 41, and 
4. 65 by 5. 15. The largest scintillator contained NE 213, and the other two contained 
NE 218 scintillator solutions. 

The scintillators were mounted on RCA 8575 photomultiplier tubes. A modified 
Owens type pulse shape circuit as described in reference 5 was used. This circuit is 
used with a two-parameter pulse height analyzer to separate proton recoil and y- ray 
events occurring in the scintillator. The proton recoil region of the two-parameter 
data set is summed by using the PREJUD computer code (ref. 6). Commercially avail- 
able electronics were used. A circuit diagram of the spectrometer system used is 
shown in figure 1. The proton recoil spectra were measured from 200 keV to a little 
above 10 MeV proton energy by using four amplifier gain settings to provide adequate 
overlap and spectrum detail. Proton recoils due to room scattered neutrons were ac- 
counted for by using a water shadow cone. 

pp 

All three scintillators were calibrated with a sodium-22 (Na ) gamma ray source. 

The energy used for the half height of the Compton edge was 1. 12 MeV. A calibration 

curve for the 1. 22-centimeter-diameter scintillator is shown in figure 2. Also shown 

fiO 

in the figure is a cobalt-60 (Co ) distribution for this scintillator which shows that the 
fiO 

two Co Compton edges are resolved. In figure 3 a low-energy calibration curve is 
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shown for the 1. 22-centimeter-diameter scintillator. The americium-241 (Am ) 
source used has a y- ray at low energy so that a photopeak is obtained by the organic 
scintillator. The resolution obtained from this photopeak at 0. 060 MeV is 34 percent. 

Since the spectrometer is calibrated with Compton recoil electrons provided by a 
convenient y-ray source and then used to measure protons, the relation between the 
pulse heights of these two particles must be known. The relation used in the present 
work was measured by Verbinski, Burrus, Textor, Love, Zobel, and Hill (ref. 2); 
however, data were obtained in the present work at 2. 8 -MeV proton energy and the re- 
sults are compared with the results of reference 2 in the section RESULTS AND DIS- 
CUSSION. These measurements were made by using the D(d, n) neutron reaction as a 
source of monoenergetic neutrons. The spectrometer was placed 1 meter from the tar- 
get of a 150-kilovolt accelerator, and the proton recoil spectrum due to the 2. 8-MeV 
neutrons was measured with the 1. 22-centimeter-diameter by 1. 27-centimeter-long 
scintillator. Background due to room scattered neutrons was accounted for by using a 
shadow cone. These data are shown in figure 4. The proton recoil distribution is 
roughly rectangular in shape as expected but rises with decreasing energy because of 



the nonlinearity of proton pulse height as a function of energy. The shape would be rec- 
tangular for an ideal spectrometer. 

DIFFERENTIAL SPECTRAL ANALYSIS 

For an ideal scintillator that is thin enough so that multiple scattering is unimpor- 
tant and yet thick enough so that wall effects are not important, the measured proton 
recoil distribution is a relatively simple function of the neutron spectrum and can be 
inverted to obtain the spectrum. Small corrections can then be applied to account for 
the deviation of the actual scintillator from an ideal one. 

The differentiation method is based on neutron-proton isotropic scattering for which 
all proton recoil energies less than or equal to the incident neutron energy are equally 
probable. For monoenergetic neutrons the number of recoil protons per unit energy per 
second dN /dE is given by 

r r 

dN 

— 2 - Ae (1) 

dE E 

P n 

where cp is the number of neutrons per square centimeter per second at energy E 

2 n 

(MeV) incident on a scintillator of area A (cm ) normal to the beam, and e is the in- 
cident neutron scattering probability of the hydrogen in the scintillator or efficiency of 
the scintillator. When polyenergetic neutrons are incident on the scintillator, the recoil- 
proton pulse height distribution dN /dE has contributions from all neutrons greater 

r . r*' 

than or equal to Ep. Therefore, 


dN 

— e = a 

dE p 


/ OO 

e(E n ) 

-n= E p 


V(E n ) 


dE. 


E_ 


n 


( 2 ) 


The neutron spectrum is obtained by differentiating equation (2) with respect to Ep and 
solving for <p(E n ): 


(pi E n ) = - 


E , / d 2 w p\ 

dE? / 

' P/ E =E 


( 3 ) 
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Since the scintillator light output does not vary linearly with recoil proton energy, the 
quantity measured is not dNp/dEp but dN^/dEg, where Eg is the energy of beta 
particles which are used to calibrate the spectrometer. Equation (3) must then be re- 
written: 


n 


<p(EJ = - 

Ae(Ej 


dN dE 


n' 



B 


dE„ dE 
B p 


Je =E 
p n 


( 4 ) 


The variation of E„ with E-r> used in this work was measured in reference 2. The 
P a 

measurements in reference 7 indicate that NE 213 and NE 218 are essentially the same 
with respect to the variation of Ep with Eg. 


Multiple Scattering and Wall Effect Corrections 


In a real scintillator there are multiple neutron scatterings, wall effects, and carbon 
interactions. 

We have accounted for multiple scattering and wall effects by using the analysis of 
reference 8 in which two corrections are derived, L for proton leakage or wall effect 
and S for multiple scattering. 

The leakage plus multiple scattering correction term is computed from the formula 


L + S = 


1 - 0.78 


R(EJ 


n' 


+ 0. 090 N h rO<E n ) + 0. 077 N h rcr(0. 68 E n ) 


(5) 


2 

where R(E_) (mg/cm ) is the range of a proton which receives the full neutron energy 

n 2 24 3 

E n (MeV), r (mg/cm ) is the scintillator thickness, (atoms X10 /cm ) is the hydro- 
gen atom density in the scintillator, cr(E n ) (barns) is the neutron scattering cross sec- 
tion for hydrogen, and r (cm) is the radius of the scintillator. The first term of equa- 
tion (5) accounts for the wall effect on the basis of 0. 78 of that part of the scintillator vol- 
ume in which protons can strike the wall. The correction is used as a flux divisor. 

The calculated correction (L + S) is plotted in figure 5 for the 1. 22- and 4. 65- 
centimeter-diameter scintillators. For the small scintillator, the correction factor 
varies from 1. 07 at 0. 5 MeV to 0. 91 at 10 MeV. The correction (0. 07), therefore, is 
small enough to allow the method of differential spectral analysis to be applied without 
resulting in a significant error in the measured spectrum. For the large scintillator, 
multiple scattering dominates at all energies, and the correction factor is large, reach- 
ing 1. 26 at 0. 5 MeV. Unfortunately, the multiple scattering correction is probably the 
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least accurate of the two corrections because the nonlinearity of the proton light output 
causes these summed pulses from multiple scattering to appear significantly below the 
leading edge of the recoil distribution where the correction is applied. Therefore, the 
large scintillator requires calibrated monoenergetic neutron response functions to 
measure spectra accurately. 


Carbon Interactions 

Carbon interactions are the most difficult to account for in this type of analysis be- 
cause of their effect on scintillator efficiency, as discussed in this section. However, 
pulses produced in the scintillator by carbon recoils and alpha particles from the C(n, a) 
reactions can be identified by pulse shape analysis and eliminated from the data. In the 
case of fission-like neutron spectra, the number of carbon reactions is small, and re- 
coil pulses fall below most lower bias settings. Here the number of alpha particle 
pulses is small and can be ignored because of the high threshold for the C(n, a) reac- 
tion (6 MeV). 

The characteristic shape of the proton recoil distribution can also be affected by 
neutrons interacting with carbon because a neutron can lose up to 28 percent of its 
energy in a carbon collision and then produce a recoil proton. The net effect of these 
processes and the hydrogen multiple scattering processes appears as an inflection in the 
proton recoil distribution at about 80 percent of the maximum energy. In the case of the 
data in figure 4, this inflection is not very evident within the statistics of the measure- 
ment, but it is very evident in the distributions shown in reference 2 for their larger 
scintillator. Since carbon multiple scattering processes, in which the collision involves 
only a small energy transfer, are indistinguishable from single hydrogen collision neu- 
trons, they should be included in the scintillator efficiency. Proton recoils following 
large energy transfer carbon collisions cannot fall in the leading edge of the distribution 
and should not be included in the efficiency calculation. 


Scintillator Efficiency Calculation 

The efficiency of the scintillator e is the probability that an incident neutron enter- 
ing the detector will create a recoil proton. Since there is some self- shielding in all 
scintillators of practical size, the efficiency will vary depending on the angular distribu- 
tion of the incident neutron flux. 

Exact expressions for e can be obtained for special cases. For the case of a par- 
allel beam of neutrons incident on the flat face of a scintillator of thickness t and area 
A, the total number of first collisions that are scattering collisions with hydrogen is 
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Ae = All - e 


( 6 ) 




where S T is the total macroscopic cross section of the organic liquid and Ejj is the 
hydrogen macroscopic cross section of the organic liquid. For the case of a parallel 
beam of neutrons incident on the curved surface of a cylindrical detector, a similar ex- 
pression exists that involves an integral over the variable thickness. 

However, neutron interactions with carbon alter the collision probability with hy- 
drogen. The scintillator efficiency including carbon scattering contributions to the pro- 
ton recoil distribution can be accurately computed from the following expression from 
reference 9: 


Ae = All - e 


-2. t 
tr 


'H 


S H + S ctr^o 


(7) 


where E_ tr is the transport cross section, 2 ctr is the carbon transport cross section, 
and p Q is the flat source escape probability (ref. 10). In expression (7), it is assumed 
that all carbon scattered sources are uniformly distributed. The expression was orig- 
inally derived in reference 9 for sphere transmission analysis. In equation (7), hydro- 
gen scattering by carbon atoms is treated as an absorption since further multiple scat- 
tering occurs within the time resolution of the spectrometer and simply increases the 
pulse size. An alternative to using equation (7) to calculate the total efficiency is 


Ae 



( 8 ) 


which completely ignores the presence of carbon in the scintillator. This expression is 
indicated in reference 1 to give accurate results for scintillators up to 1 centimeter 
thick. 

For the case of the scintillator immersed in an isotropic flux, the efficiency can be 
calculated by using methods developed for calculating resonance neutron capture in 
thermal reactors. Here the average fluxes in the two adjacent media are obtained by 
solving the coupled integral equations for the energy dependent fluxes that express the 
neutron balance in each region. The spatial transfer of neutrons between the adjacent 
media are obtained by means of their respective collision-escape probabilities p Q and 
application of a reciprocity condition. A large amount of experimental and analytical 
information is available on approximations to p Q and the accuracy of this method 
(ref. 10). The expression for the proton production rate per unit flux Ae is 
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(9) 


Ae = p Q S H At 


for the first collision efficiency or 


Ae=p 0 S T At £ (10) 

S H + p o S ctr 

for the total collision efficiency, where At is the volume of the scintillator. 

Some insight into the applicability of various sized scintillators to differentiation 
analysis can be obtained from the previous expressions. Figures 6 and 7 show calculated 
efficiencies for various conditions as described in the figures. 

In figure 6, the calculated efficiency times area is shown for a 4-65-centimeter- 
diameter by 5. 11-centimeter-long NE 213 scintillator. The upper curve is for the case 
of a parallel beam incident on the side of the scintillator and includes multiple scatter- 
ing contributions. The three lower curves are for the types of flux incidence shown but 
do not include multiple scattering. Inclusion of multiple scattering increases the effi- 
ciency of this size scintillator approximately 15 percent at 1 MeV. As mentioned 
earlier, because proton recoils from carbon collisions cannot fall on the leading edge of 
the pulse height distribution, it is not clear that all of this calculated increased efficiency 
is applicable in differentiation analysis. Therefore, the true efficiency is bounded by 
the two upper curves. At higher energies this variation reduces to only about 5 percent. 
This energy dependent aspect of the calculated efficiency will give a small shape error 
to the derived spectrum. The calculated efficiencies in figure 6 also show a difference 
at 1 MeV of 18 percent for a parallel neutron beam incident on the flat surface compared 
to one incident on the curved surface of the scintillator, so that measurements with the 
larger scintillator are indicated to be sensitive to the angular distribution of the incident 
flux. 

Similar efficiency calculations are shown in figure 7 for a 1. 22-centimeter- 
diameter by 1. 27-centimeter-long scintillator. The calculations show that multiple 
scattering is significant even for this small scintillator but that the efficiency variation 
involved is less than 5 percent at 1 MeV. The sensitivity of the efficiency of the smaller 
scintillator to the angular distribution of the incident neutron flux is also less than 
5 percent for neutron energies greater than 1 MeV, which shows that flux measurements 
can be accurately made where the angular distribution is unknown. The efficiency is _ _ 
also low enough so that measurements can be made in a reactor. 

Equations (4), (5), and (8) in the preceeding analysis are incorporated in the com- 
puter code TUNS written for the IBM 7094 computer. Parallel-beam neutrons incident 
on a curved surface can also be treated. A description of the code is given in appen- 
dix A, and the program listing is given in appendix B. 
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RESULTS AND DISCUSSION 


Prior to comparing spectrum measurements of the PoBe source, it is of interest 
to analyze the proton recoil distribution for the 2. 84-MeV monoenergetic neutron source 
presented in figure 4. Generally, the Compton electron energy corresponding to the 
half height of the leading edge of the proton recoil distribution is used as the beta energy 
corresponding to the neutron source energy; from the data in figure 4, this is 0. 96 MeV. 
The calibration data for Ne-213 in reference 2 indicate that a 0. 96-MeV beta corre- 
sponds to a proton energy of 2. 78 MeV; this is 2 percent lower than the known neutron 
source energy of 2. 84 MeV but within the 2 percent error quoted in reference 2. 

There is some uncertainty in determination of the half height of the data in figure 4 
since the full height is only characterized by a marked reduction in slope. An alterna- 
tive method of determining the peak flux is to apply differentiation analysis to the data. 
The results show that the peak in the spectra can be fairly well determined but that very 
narrow slope-taking intervals must be used in the vicinity of the peak. It is therefore 
important that good counting statistics be obtained in order to use this method. The 
neutron spectra also show that between 2. 0 and 2. 4 MeV there is an indication of low 
residual flux values which may be attributed to multiple scattering in the scintillator. 

The resolution of the spectrometer can be obtained from figure 8 and is 8. 6 percent 

at the source energy of 2. 84 MeV. The resolution at this energy of the 5-centimeter- 

diameter by 5-centimeter-long scintillator used at ORNL is given in reference 11 and 

is 20 percent. A second comparison can be made by using the data Of figure 3, where 
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the Am photopeak data which correspond to a neutron energy of 0. 47 MeV show a 
resolution of 34 percent for the 1. 22-centimeter-diameter scintillator. The resolution 
at this energy given in reference 12 for the 5- by 5-centimeter scintillator is 50 percent. 
The poor resolution reported shown in reference 12 is due in part to the inclusion of 
adequate smoothing in the FERDOR analysis of the measured data. 

A significantly better resolution is indicated for the smaller scintillator. Although 
more data are desirable, the two-point comparison indicates a difference in shape of 
resolution as a function of energy for the two spectrometers. 

The measured proton recoil spectra for the PoBe source using three scintillators 
are shown in figure 9. The small scale of the figure tends to smooth out much of the 
detail actually in the data, but some of the larger differences in the shapes of the three 
distributions are evident. Below a proton energy of 3 MeV (1 MeV beta) the curve for 
the 1. 22-centimeter scintillator rises more rapidly than that of the larger scintillators 
because of less neutron self-shielding. At a proton energy of 10 MeV (5 MeV beta), a 
more pronounced inflection in the curve for the small scintillator is apparent because of 
its superior resolution. 

Figures 10 to 12 show the results of the differentiation analysis for the PoBe source 
measurements. Figure 10 shows the spectrum obtained by using the NE 213 scintillator 
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that was 4. 65 centimeters in diameter and 5. 15 centimeters long. The data were ana- 
lyzed by using the differentiation method and the unfolding method with the FERDOR code 
(ref. 13) and the data of reference 2. The FERDOR results were reduced 12 percent to 
correct for the greater length of our scintillator and to make our results directly com- 
patible with the scintillator size used in reference 2. The spectrum obtained by differ- 
entiation is seen to be significantly lower than the spectrum generated by the FERDOR 
code, throughout the range of energy. The same proton recoil pulse height distribution 
data were used for both codes. The shapes of the two spectra are in general the same 
in both cases, except for a difference in energy location of the 3. 4-MeV peak indicated 
by FERDOR and the presence of a peak at 2. 2 MeV shown by the FERDOR results; a 
corresponding change in slope in the measured proton recoil distribution was not present 
at 2. 2 MeV. 

The deep valley that exists in the PoBe neutron spectra at approximately 1. 5 MeV 
is a good test of proton recoil data differentiation analysis, since any deviation of the 
proton recoil distributions in the vicinity of 3. 4 MeV from ideal shapes will be apparent 
in the unfolded spectra near 1. 5 MeV. Thus, the change in slope of 3. 4-MeV distribu- 
tions at 80 percent of maximum pulse height due to carbon interactions, which are not 
accounted for by differentiation, could combine with small slope changes due to 2. 2-MeV 
neutrons and result in a single slope change to produce the peak in the data at 3 MeV. 

Figure 11 shows the spectrum obtained with the 2. 16 -centimeter-diameter scintil- 
lator. Also shown in the figure for comparison are the FERDOR results for the 4. 65- 
centimeter-diameter scintillator. The uncertainty for the FERDOR results is not in- 
cluded for clarity. The differentiated results for this size scintillator show much better 
agreement with the FERDOR results than with the differentiated results for the 4. 65- 
centimeter scintillator; however, there is still not a good match in the shape of the 
spectrum in the vicinity of 2 MeV. 

Figure 12 shows the spectrum obtained with the 1. 22-centimeter-diameter scintil- 
lator; this spectrum is also compared with the absolute FERDOR spectrum shown in 
figure 10. Comparison of the two spectra indicates good agreement with clear indica- 
tions of the better resolution of the 1. 22-centimeter-diameter scintillator. The small 
peak at 2 MeV indicated by FERDOR is also suggested in the data of figure 12. A second 
method of intercomparing the spectra is to compare the absolute integral flux. This has 
the advantage that the integral should be the same for two measured spectra even though 
the respective resolutions are different. Table I shows the integral flux above 0. 5 MeV 
for the four determinations. These results show that the FERDOR integral flux is the 
same as that obtained for the two smaller scintillators but 15 percent higher than that 
obtained by differentiation for the 4. 65-centimeter-diameter scintillator. The statistical 
uncertainty obtained with FERDOR is about one -half that obtained by differentiation. The 
statistical uncertainty is higher for the smallest scintillator because of insufficient 
counting times used in obtaining the proton recoil distribution. 
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The calculated total efficiency for proton production was used in the differentiation 
analysis. This efficiency includes protons produced following carbon collisions; there- 
fore, it is not too surprising that the resulting flux for the 4. 65-centimeter-diameter 
scintillator is lower than the small scintillator and FERDOR results. This comparison 
gives some credence to the postulate that only a small fraction of the protons produced 
following carbon collisions should be included in the scintillator efficiency used in the 
differentiation analysis, and that an effective efficiency for use in this analysis is the 
first collision efficiency. Unfortunately, part of the essentially constant percentage dif- 
ference between the magnitude of the flux can also be explained by a large error in the 
L + S correction to the analysis, which also would result in low fluxes, as in the sec- 
tion Multiple Scattering and Wall Effect Corrections. 


IMPLICATIONS OF RESULTS 

Measurements have been made of the PoBe neutron spectrum by using three organic 
liquid scintillators of different sizes. Measurements were also made for 2. 8-MeV neu- 
trons from a D(dn) source by using a small 1. 22-centimeter-diameter by 1. 27- 
centimeter -long scintillator. An inter comparison of the absolute spectra was made 
from a calibrated unfolding analysis for the largest scintillator and a differentiation 
analysis for all the scintillators. Results indicated that spectra consistent in shape and 
magnitude can be obtained by differentiation for small liquid scintillators between 1 and 
2 centimeters thick. If differentiation is used for scintillators as thick as 5 centime- 
ters, significant shape differences in spectral valleys can result. The proper effective 
scintillator efficiency to be used in the analysis is indicated to be close to the first 
collision efficiency. 

In the case of the smallest scintillator, the efficiency is shown to be insensitive to 
multiple scattering and the angular distribution of the incident flux. This characteristic 
and the inherently superior resolution of this size scintillator make its use desirable 
for flux measurements within diffusion media and in applications where high efficiency 
is not required. A computer code TUNS was written to analyze the data by differentia- 
tion analyses. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 20, 1972, 

503-10. 
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APPENDIX A 


COMPUTER CODE TUNS 
Code Description 

The computer code TUNS which solves numerically equation (4) has been written in 
FORTRAN IV for the Lewis IBM 7094-n computer. The value of a derivative at a point 
is obtained as the average of the slopes from the point to its nearest neighbors, weighted 
inversely as the distance between points. Analytic expressions are used to compute 
R(E n ), a (E n ), and a(0. 68 E n ) in equation (5) (ref. 14). The efficiency is calculated by 
using equation (8). A function subroutine SIMPS 1, (ref. 15) is used for integration in 
computing e for a parallel beam of neutrons incident on the curved surface of the scin- 
tillator. 

Input is described in the next section. The program output is in terms of neutron 
energy, flux, and the statistical error. Errors are computed from the standard devia- 
tions applied as input. The energy spacing between points is chosen to correspond to 
the spectrometer resolution. A wide spacing can be used and a small error will result, 
but spectral detail will be lost. The effect of the size of the interval from which the 
slope is determined bn the accuracy of the spectra is discussed by Bennet (ref. 16). 


Input Format for TUNS Code 


Card 

number 

Card 

column 

Format 

Variable 

name 

Description 

1-14 

1-80 

8G10. 8 

PE 

dimensioned variable; should con- 
tain Ep in increasing order 

15-28 

1-80 

8G10. 8 

DLDEP 

dimensioned variable; should con- 
tain the values of dE^/dE,, 

a p 

evaluated at Ep points for PE 
in increasing order of Ep 

29-42 

1-80 

8G10. 8 

LT 

dimensioned variable; should con- 





tain light output in terms of 
cobalts 

43 

2-73 

12AG 

TITLE 

title card describing experiment 

44 

1-12 

F12. 5 

DEBDL 

conversion factor used in this 
report, 1.240 MeV/Co 
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Card 

Card 

Format 

Variable 

Description 

number 

column 


name 



13-24 

F12. 5 

FACE 

1. 0 if flux is incident on flat sur- 
face 


28-30 

16 

NPTS 

number of input values (see card 4) 

45 

1-12 

F12. 5 

NH 

hydrogen atom density, 
atoms xlO^/cm^ 


13-24 

F12. 5 

R 

radius of scintillator, cm 


25-36 

F12. 5 

H 

height of scintillator, cm 


37-48 

F12. 5 

BW 

analyzed channel width, MeV 


49-60 

F12. 5 

ADJ 

analyzer live time, sec 


61-72 

F12. 5 

WREF 

1. 0 if A, L, and S are desired as 
output 

46(1) 

26-28 

15 

N(I) 

card number 

1=1, . , NPTS 

41-60 

E20.9 

DNDC(I) 

counts per channel 


61-80 

E20.9 

ERR(I) 

standard deviation in counts per 
channel 


Card 46 (NPTS) may be followed by another set of data, beginning with a card 43; as 
many sets as desired may be run. Cards 46(1) to 46 (NPTS) may be obtained from 
code PREJUD (ref. 10). 
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APPENDIX B 


PROGRAM LISTING 


STOP TIM£ = 1 ,PAGES=15,FURMS=PLAIN 

S I 8 JOB DEBUG 

SlBFTC TUNSBC DEBUG 

DIMENSION PEI 109 I .EN( 1091 .OLDEP 1 1C 91 . DLOE (109 ) ,ONDC 1 109 I « ARG( 1C9 ) » 
1UARGI 10 71 .FLUX! 107 ) , N 1 1C9 ) .ERR (109 ) .RELERR 1 109 I ,F LUP ( 107 » • AER ( 1*39 ) 
2 « FLO ( 107) , SFLUX ( 1C 9) » UF L U X ( 109) 

RFAL LT (109). LG( 1C9) .NH 
COMMON S.SIG 
EXTERNAL TRANS 
READ! 5 » 1000 ) PE 
R E AC ( 5 • 1000 ) ULDEP 
READ! 5. 1000) LT 
l R EAC ( 5 » ?50Q ) 

WRI T E I 6 ,2400 ) 

WR l T E ( 6 ,2 500 ) 

R E AC ( 5 « 1 1 5C ) DEBDL, FACE, NPTS 
READ! 5. 110C)NH*R*H*8w*ADJ » WREF 
READ! 5, 1200) (Nil ),DNDC( I I, ERR ( I I ,1 = 1, NPTS) 

WRITE (6,10 50 I OEBDL , NH ,R«H*&W»ADJ 
WR I T E ( 6 , 1060 ) NP1S 
I F( FACE-1. O) 3,4,3 

3 WRITF(£.1C75) 

REF =S0R T ( 2.0*R*H/3.1416I 
GO TO 5 

4 WRITE (6 ,1080 ) 

R EF=R 

5 WRITE(6,1085) 

WRITE (6 ,1090) ( t,N( I) ,DNDC(I) ,ERR( I) , 1=1, NPTS) 

00 10 1*1, NP T S 

(F(N( I ) ,E0.28.0R.N( I).E0.30.0R.N( I >.E0.32.0R.N( I) .E0.34) GG TO 6 
GO TO 7 

6 WR I TF ( fc ,1070 ) N ( I ) 

STOP 

7 I F ( N ( I) .E0.29) N( I)=28 

1 F( N ( ll.E0.31) N ( 1 ) =29 
(F(N< II .E0.33 I Mil =3C 
IF(NI I ) .GT.34) N ( 1 1 *N( 1 1—4 
IF(N( I) .GT.1G9I GO TO 8 

GO TO 9 

fl WRI TE ( 6 ,1065 ) 

STOP 

9 M=N( I ) 

FN( I )=PE(MI 
OLOE( I )*OLDEP(M) 

LG(l)=LT(M) 

10 CONTINUE 
12 00 14 I ® 1 » NPT S 

ONf)L = OEeOL*ONOC( I )/BW 
ARG( I )=CNDL*DLDE( I I 
AER ( I ) = CL DEI I ) WOE BDL*ERR ( D/BW 
14 CCNTINUF 

MPT S=NPTS-2 
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15 DO 20 1 =1 • MPTS 

0 A= ( AR C- ( I +1 )— ARG I 1 ) ) / ( LG ( i+l)-LGt I ) ) 

0B=(AKG(I+2)-AKU( I + ll )/(LG( l+2)-LG< 1+1) ) 

DmA =( L C- < I +2 I -LG ( 1 + 1) )*CA 

Ot«fi = ( L C- 1 l + 1 )-LG I 1 II *06 

OARGI l ) = < OwA+DwB ) /( LGl I *2 l-LG ( I ) ) 

RELERRI n=A8S( SCRT I AER ( I ) **2+AER< 1+2 ) **2 ) / ( ARG( I ) -ARG ( 1+2) ) ) 

20 CONTINUE 

IFIMREF-l.O) 2 5.24.25 
2A MR l T E ( 6 . 1600 ) 

2 5 OG 50 [ =1 . MPT S 

HOT = 1.0+7.A17*EN( l+l ) +C.1 1C5*EM l+l)**2 
E IR ST = 1 7. 602 A/BCT 

B0TT=1 .C+0.2A27*ENI I + 1 ) +0 .002 6*EN (1 + 1 ) * *2 
SEC-2.7181 /80TT 
S IGH = F IRST+SEC 
S I G=NH* S I GH 
X =0 

IFIFACE-1.0) 27,25.27 
27 SlM=SIMFSl(O.C.R. TRANS. Kl 
FF = R-S IP 
A=2.0*R*H 
FFF=EF/R*A 
S IMR = S I P /R 
Th=ALOG IS I MR) 

T HFF =— 1 .0*TH/SIG 
GO TO 2 1 

29 FF=1 .C-EXPI-1.C*SIG*H) 

A=3 . 1 A 1 6* H*H 
EFF=EF* A 
TFFF=H 

31 ENL=G.C68*EN< I +1 ) 

D0MN = 1 .C+7.A17*ENL+0.11C 5*ENL **2 
0NE=1 7.6G2A/OUMN 

OCMNN= 1 .0 + 0. 2A27*EML+C.C028*ENL**2 
T i»U = 2 . 7181 /DOWNN 
S IGHL = CI\E + TWU 

RM=1 .73F2*(ENI I +1 ) +C . 1 50 A 5 ) ** 1. 81 9 A 
RANGE-RP/870. 

ESC = 1 .0-0. 78Q*RANGE/TFEF 
SESA=0.090*NH*THEF*SIGH 
S ES 8 = 0. < 7 7 *NH*RF F*S l GHL 
S£SC=SFSA+SESB 


8EE=FSC +SESC: 

IF) WRFF-l.C) 36.2A, 36 

3 A wRI TEI6 .1700) EM i + 1 ) .EFF ,K .ESC, SESC , RANGE 

36 0N0F=-1 .0*EN( 1+1 )*QLDt< I+1)*0ARG( I ) /EFF /BEE 
FLUX ( I ) =GNOE/AGJ 
PREG = R E LER R ( I )*FLUX ( I ) 

ELUPI I )=FUJX< ( ) + PREC 
Fl.O( I ) = FL U X I I l-PREC 

50 CONTINUE 

IM = 0 

SFI.UXI 1 ) = 0. 

Ilf LUX < 1C8) = G. 

no 70 I = 2. MPTS 
.) = MPTS + 2-1 

SFLUX(I) = SFLUXII-1) ♦ FLUX(I) * (ENII+1) - EN(I)) 
UFLUXIJ) = UFLUXIJ+1) + FLUX(J) * (EMJ+1) - EM J ) ) 



IF (FLUX! I) .GE. 0.) GO TO 70 

IM = I 

SFLUX! I ) = 0. 

70 CONTINUE 

IF ( IM .EG. 0) GO TO IOC 

DO 80 I = 1. IM 

SFLUXI 11=0. 

80 UFLUXI I » = UFLUXI IM + 1 ) 

IM = IM + 1 

SFLUXI IM = FLUX I I M ) * IcKlIIM + 1) - EN ! I M) ) 

100 WRITEI6.1500) 

WRITE (6,20001 ( UNI I+lt.EM Ul) • FLUX ( 1) ,FLUP( I ) ,FLO( I ) , SFLUXI I), 

♦ UFLUXI I I. I = l.MPTS) 

GO TO 1 

1000 FORMAT (8G10.8I 

1050 FORMAT! 1H0 , 66HFACT0R FOR CONVERTING COUNTS PER MEV BETA TO COUNTS 
1 PER COeALT IS , 1PE12. 5/ /3GH THE HYDROGEN ATOM DENSITY IS ,1PE12.5/ 
2/35H THE RADIUS OF THE SCINTILLATOR IS .1PE12.5//35H THE HEIGHT OF 
3 THE SCINTILLATOR IS . 1PE 12. 5 / /28H THE BIN WIDTH (M VALUE! IS ,1PE 
A12.5//31H THE FLUX ADJISTMENT FACTOR IS .1PE12.5! 

1060 FORMAT! 1H0.30HTHE NUMBER OF INPUT VALUES IS ,13) 

1065 FORMAT! 1H0 , 35HERROR- LOC. NO. IS GREATER THAN 109) 

1070 FORMAT! IHC.52HERR0R- CARD SHOULD BE REMOVED FROM DATA. CARD NO. = 
1,13) 

1075 FORMAT! 1H0.62HTHE FLUX IS INCIDENT ON THE CURVED SURFACE OF THE SC 
1 INT ILLATOR ) 

1080 FORMAT! 1H0.60HTHE FLUX IS INCIDENT ON THE FLAT SURFACE OF THE SCIN 
IT ILL ATOM ) 

1085 FORMAT! IHO ,6X , 9H INPUT NO. ,10X , 8HCAR0 NO. , 20 X , 6H COUNT S*1AX,9HSTD. D 
1EV.) 

1090 FORMAT! EX. 13 . 18X . 1 3 , 20 X , 1 PE 12 . 5 . 8X . 1PE1 2. 5 ) 

1100 FORMAT (6F12.5) 

1150 FORMAT! 2F12.5. 16) 

1200 FORMAT! 23X, 1 5 . 12X ,2E20. 9 ) 

1500 FORMAT ( 1HL .BOX , 2 1 10H l NTEGKAT ED . AX) / AX , 7HPT. NO. , 2X, 8HL0C . N0..8X, 
16HENERGY. 10X, 4HFLUX ,?X. llHUPPER BOUND ,3X. 11HLUWER BOUND, 2X , 12HFLUX 
2 BELOW E« 2X , 12HFLUX ABOVE E// ) 

1600 FORMAT! }H0 , 6X . 6HENERG Y , 13 X , 10FEFF IC I ENC Y , 6X .4HERRP ,21 X , 7HLEAK AGE , 
112X.7HSFCSCAT «12X« 5H RANGE ) 

1700 FORMAT! ? ( AX , 1PE 12 . 5 ,4X ) , I 5,1 5 X ,3! AX , 1PE 12 . 5 , 4X ) ) 

2000 FORMAT! 8X, 13, 1 10 , 1P6E1A. 5 ) 

2A0C FORMAT ( 1H1 ) 

2500 FORMAT! EOH 
1 

END 

1.IBFTC TRANSM 

FUNCTION TRANS(X) , 

COMMON P.SIG 

PL=-2.0*SIG*SUKT!R**2-X**2> 

T RAN S = E XP ( PL ) 

RETURN 

END- ----- * - 
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TABLE I. - MEASURED TOTAL FLUX ABOVE 


0. 5 MeV AT 1 METER FROM POLONIUM- 
BERYLLIUM SOURCE 


Scintillator 

size, 

cm 

Unfolding 

method 

Total flux, 

2 

neutrons/(cm )(sec) 

4.65 by 5. 15 

FERDOR 

99. 0±0. 8 


Differentiation 1 

85. 8±1. 7 

2. 16 by 2.41 | 

Differentiation 

99. 1±1. 9 

1.22 by 1.27 

Differentiation 

103±2. 8 


RCA8575 



JO02m^ 
High 
voltage 


10.01 mF 


Figure 1. - Spectrometer block diagram and photomultiplier tube base circuit. 










Figure 3. - Americium-241 calibration curve for 1. 22-centimeter- 
diameter scintillator. 


.8 


.9 


Energy, MeV 

Proton recoil distribution from 2. 8-MeV neutrons using 1. 22-centimeter-diameter scintillator. 


Scintillator dimensions, 
cm 


4. 65 by 5.11 
1. 22 by 1. 27 
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